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We investigate the effect of (weak) dipolar interactions on the field behavior of the 
temperature at the maximum of the zero-field-cooled magnetization of a polydis- 
perse assembly of nanoparticles. For this purpose we extend the Gittleman-Abeles- 
Bozowski model for the zero-field-cooled magnetization by computing the contribu- 
tion of dipolar interactions to the longitudinal relaxation time. We show, in good 
qualitative agreement with many experimental observations, that the temperature 
at the maximum of the zero-field-cooled magnetization as a function of the applied 
field changes from a bell-like to a monotonically decreasing curve when the intensity 
of the dipolar interactions, or equivalently the sample concentration, increases. 

PACS numbers: 75.50.Tt; 75.10.Hk; 05.20.-y 



I. INTRODUCTION 

One of the direct technological applications of magnetic nanoparticles is magnetic record- 
ing. The storage capacity of the media can be considerably increased by devising denser 
assemblies of smaller and smaller particles. However, this brings a dilemma because small 
particles become superparamagnetic, i.e., thermally unstable, well below the room temper- 
ature. Moreover, high density in assemblies entails strong dipole-dipole interactions (DDI) 
between the particles, and in technological applications such as magnetic recording, this is 
an issue of special importance since DDI have been widely recognized as being responsible 
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for the deterioration of the signal-to-noise ratio (see e.g., Refs. Ill, |2| and references therein). 
As such, an optimal material, with appropriate anisotropy and other physical parameters, 
has still to be devised. On the other hand, the study of nanoparticle assemblies brings new 
hurdles to theorists, at least, since one is faced with tremendous difficulties related with 
DDI between particles, together with the distributions of volume and anisotropy axes. In 
spite of that, DDI in nanoparticle assemblies has triggered much interest due to many new 
phenomena that emerge from the collective behavior of the particles and also because these 
interactions have always constituted a challenging issue in many areas of physics. 

To analyze and eventually understand experiments on the dynamics of interacting 
nanoparticle assemblies, and in particular to understand the dynamical response, such as 
the ac susceptibility and zero-field-cooled (ZFC) and field-cooled (FC) magnetizations, one 
needs to know how DDI affect the switching process and the relaxation time of the nanopar- 
ticles. While a fair understanding of the mechanisms underlying the ZFC-FC magnetization 
process and ac susceptibility has been achieved in the case of non-interacting assemblies, 
many experimental results on interacting assemblies remain unexplained. Obviously, this is 
mainly due to the long range of DDI and also to the complexities of the very calculation of 
the relaxation time itself. Recently, Jonsson and Garcia-Palacios [see also {4]] obtained 
an approximate expression for the relaxation rate of a weakly interacting monodisperse 
assembly of macropins with textured or randomly-distributed anisotropy. The macrospin 
approximation here means that a particle is represented by a macroscopic magnetic moment, 
i.e., ignoring its internal structure. This picture of the particle will henceforth be referred to 
as the one-spin problem (OSP). In the literature, this is also known as the coherent-rotation 
limit. Then using the simple Debye relaxation model they investigated the effect of DDI 
on the ac susceptibility and in particular the displacement of the maximum of its real and 
imaginary components. In their explanation of this effect they emphasized the important 
role played by damping in the relaxation processes in the presence of a transverse field in 
addition to the effect of the change in the energy barriers, that was commonly believed to 
play the major role. The role of a transverse field is played here by the transverse component 
of the dipolar field. 

Experimental results obtained for ferrofluids 5| and later for 7~Fe203 nanoparticles 6] 
indicated that for dilute samples (weak DDI), the temperature T max at the maximum of 
the ZFC magnetization first increases with increasing field, attains a maximum and then 
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decreases. More experiments performed on the 7-Fe 2 3 particles dispersed in a polymer 0,0] 
matrix confirmed the previous results for dilute samples and showed that, on the contrary, 
for concentrated samples (strong DDI) T max is a monotonically decreasing function of the 
magnetic field. The shift of this maximum was also studied with different techniques in 
various types of nanoparticles, see for instance j^], Q, 11]. In Ref. Isl it was shown that the 



bell-like shape of T max (if) is not very sensitive to the intrinsic properties of the particles, of 
course in the OSP approximation. Exact numerical calculations 3, H, Q] of the smallest 
eigenvalue of the Fokker-Planck operator invariably led to a monotonic decrease in the 
blocking temperature, and thereby in the temperature T max , as a function of the magnetic 
field. Indeed, it was shown that the expression of the single-particle relaxation time does 
not play a crucial role and that even the (relatively) simple Neel-Brown expression for the 
relaxation time in a longitudinal field leads to a maximum in T max (if). What seemed to play 
a crucial role is the fact that the magnetization, formulated within the Gittleman, Abeles, 
and Bozowski (GAB) model 15], has a superparamagnetic contribution that is a non-linear 
function (such as Langevin's) of the magnetic field. The magneto-crystalline anisotropy and 
the volume-distribution width also have strong influence. The issue of the effect of DDI 
on T max (H), namely the disappearance of the maximum when the intensity of interactions 
increases, was left open in Ref. |8|. In the present work we revisit this issue after generalizing 
the work of Ref. |3| to include the static magnetic field and magnetic- moment distribution 
(polydisperse assembly) of nanoparticles in the OSP approximation. We then investigate 
the effect of (weak) DDI on the ZFC magnetization and in particular on T maiX (H). 

The present work is organized as follows: After fixing the notation and defining the 
model Hamiltonian, we explain our formalism for computing the ZFC magnetization within 
the GAB model: we compute the contribution of DDI to the longitudinal relaxation rate, 
and then show how the GAB model is accordingly extended. Next, we discuss the effect 
of DDI on T max (H) of a given assembly and compare the results for two materials, namely 
maghemite and cobalt. 



II. NOTATION AND BASIC FORMULAS 

We consider an assembly of magnetic moments rrij = mjSj, % = 1, . . . , ]\[ of magnitude rrii 
and direction s,, with |sj| = 1. The magnitude of the magnetic moment rrij is given in terms 
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of the Bohr magneton fj, B , i.e., m 8 = n^Bi and the numbers rii are log-normal distributed. 
Each magnetic moment is assigned a uniaxial easy axis e i? and for the assembly these axes 
are randomly distributed. The energy of a magnetic moment nij with uniaxial anisotropy 
axis ej, interacting with all the others via DDI, in the magnetic field H = Heh, reads [after 
multiplying by — j3 = —l/ksT], 

Si = Xi (sj • e h ) + <Ti(si • e;) 2 + f d J^ra^-Si • ^ij • Sj, (1) 

where = xrij, <7j = <7o7ij with 

_ n B H \i b K _ ///q\ //|/q 3 

being the dimensionless energy parameters. Note that <7j = KVijiksT) (or simply a for a 
monodisperse assembly) is the commonly used notation for the reduced anisotropy-energy 
barrier height of the particle %. T> is the DDI tensor defined as 

^3 = l3~ — 1) • (3) 

ij 

where = — Tj, with = r^/r^, is the vector joining the sites i and j and whose 
magnitude is measured in units of a, a characteristic distance on the matrix in which the 
particles are embedded. More precisely, the parameter a is taken as a real number times the 
mean diameter D m of the assembly, i.e., a = k x D m . Thus, large values of k correspond 
to an isotropically inflated lattice with large distances between the magnetic moments, and 
thereby weak DDI. 



III. ZERO-FIELD-COOLED MAGNETIZATION AND EFFECTS OF DIPOLAR 

INTERACTIONS 

A. Zero-field-cooled magnetization 

The dynamic response of the OSP assembly is given by the ac susceptibility. For a given 
particle with an arbitrary angle ip between its anisotropy easy axis and the field direction, 
the effective susceptibility may be written as 



XM = X|| cos2 i> + X± sin2 ^- 



Shliomis and Stepanov [16] proposed a simple Debye form for which can be generalized 
to describe the effect of a longitudinal bias field by writing 

X\\(T,H) 2 X ±(T,H) . 2 

Xshs = -r^-: cos ip + — — sin ip, (4) 

1 + iootv, 1 + zcjtj_ 

where t\\ and rj_ are appropriate longitudinal (interwell) and transverse (intrawell) relaxation 
times; x\\(T,H) and x±(T,H) are respectively the longitudinal and transverse components 
of the equilibrium susceptibility. 

In the limit of a high anisotropy-energy barrier, i.e., a ^> 1, h = x/(2a) <C 1, approximate 



expressions were found in Ref. [17| for the longitudinal and transverse components of the equi- 
librium susceptibility. If one sets in (J4j) t± = [instantaneous intrawell transverse response] 
and uses the high-barrier approximation one arrives at the GAB model [15J, generalized to 
H 7^ and an arbitrary anisotropy-axis orientation. More precisely, upon evaluating the 
high-barrier expressions for H = 0, inserting the result in (j3J), setting t_l = 0, and averaging 
over an assembly with randomly distributed anisotropy axes, one arrives at the expression 
proposed in Ref. Ila . 

Xgab ^ Xo 1 , a - -TT^ — • ( 5 ) 
The real and imaginary components then read 

, _ Xo + xA^rf) „ urq (xi - Xo) 

M?(T)V M?(T)V 
Xo 3k B T ' Xl 3KV K ] 

are respectively the susceptibility at thermodynamic equilibrium and the initial susceptibil- 
ity of particles in the blocked state (see [3] and references therein). Accordingly, starting 
from ([6]) the application of an alternating field yields: a) x' = Xo if UJT "C 1. At high tem- 
perature the magnetic moments orientate themselves on a great number of occasions during 
the time of a measurement, and thus the susceptibility is the superparamagnetic suscep- 
tibility Xo- b) x' — Xi ^ UJT ^ 1- At low temperature the energy supplied by the field is 
insufficient to reverse the magnetic moments the time of a measurement. Then, the suscep- 
tibility is the static susceptibility Xi- Between these two extrema there exists a maximum 
at the temperature T max . x' can be calculated from ([6]) using the formula for the relaxation 
time r appropriate to the anisotropy symmetry, and considering a particular volume V one 



where 
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can determine the temperature T max . Expression dSJ) of the dynamic susceptibility obtained 
for instantaneous transverse response is particularly suitable for the calculation of the ZFC 
magnetization and T max . Indeed, Eq. (JBJ) was used in Ref. si to study the effect of anisotropy 
and volume distribution on T max (if). The formalism used can be summarized as follows. 
In an assembly of particles with a volume distribution, x' can be calculated by postulating 
that at a given temperature and given measuring time, certain particles are in the super- 
paramagnetic state and that the others are in the blocked state. The susceptibility is then 
given by the sum of two contributions [jjj] 

Vc oo 

x'(T,is) = J VV X o(T,V,is) + J VV X i(T,V,u), (8) 
v c 

where T>V is the measure of the log-normal volume distribution with parameters Vq and 5 



W = — == exp 



log 



2( V_ 

V , 
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-y- 0) 



V c = V C (T, H) is the critical volume defined as the volume for which r,r = z/ m , where v m is 
the measuring frequency. V c is the "critical volume" that discriminates between the domi- 
nating populations of superparamagnetic particles of volume V < V c and blocked particles 
with V > V c , and is experiment-dependent. 

Eq. ([8]) can be rewritten for the ZFC magnetization as follows 

v c 

M zfc (H,T } ^) = Jw M sp (H,T } V,^) (10) 



o 

oo 



+ J VV M b (H,T,V^) 

v c 

where M sp = Hxo and = Hxi are the contributions to the magnetization from the 
superparamagnetic and blocked particles, respectively. 

In the present work we extend this formalism to include (weak) DDI and to investigate 
their effect on T max . For this purpose, it is necessary, in principle, to compute the contribu- 
tions of DDI to both the equilibrium susceptibility in the numerator of Eq. §5§ and to the 
relaxation time t\\ in the denominator. For the first calculation we can differentiate with 



respect to the applied field the expression obtained in Ref. [l9| for the longitudinal magne- 



tization taking account of DDI and anisotropy. We can also derive an expression for the 



transverse equilibrium susceptibility as a response to a transverse ma 



gnetic field [20] . For the 
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second calculation, we generalize the expression obtained in Ref. 
(including DDI) to include the volume distribution and the static magnetic field. This is 
done in the next section. While the outcome of the first calculation is an insignificant quan- 
titative modification of T max , the DDI contribution to the relaxation time yields a dramatic 
qualitative and quantitative effect since T max changes from a bell-like to a monotonically 
decreasing function. Indeed, with regard to Eq. ffTU]) . we will show that the change in the 
relaxation rate due to DDI induces a change in the critical volume V c and thereby a change 
in the dominating population of blocked or superparamagnetic particles. 

B. Effect of DDI on the relaxation rate 

The idea is to introduce the (local) dipolar field & (note that T>a = and see notation 
in section HIT). 

C H C/»,^>/P;;K;. (") 

3 

The relaxation rate of a magnetic moment that experiences this field can be computed using 

21l an estimation of 



perturbation theory assuming that | £j |-C 1 |3J, |4j]. Accordingly, in Ref. 
DDI was given for two samples of cobalt nanoparticles which indicates that the DDI field is 
of the order of 300 Oe, which in reduced units, obtained after dividing by the corresponding 
anisotropy field of the order of 0.3 T, is \£/H a \ ~4x 10~ 3 — 10~ 2 . This is of course very 
small, which indeed suggests that in typical (relatively dilute) samples the above condition 
on the DDI field is often satisfied. On the other hand, the magnetic field at which T max (H) 
has a maximum at approximately 100 Oe [see Fig. 1 of Ref. a, for maghemite particles], 
corresponding to a reduced field h = H/H a ~ 3 x 10~ 2 . This shows that even though £ 
is small it may still have a strong effect on T max (if) because it is of the same order as the 
applied field in the relevant range. 



The relaxation rate of a weakly DDI-interacting nanomagnet obtained in Ref. 
written as 



r 



(0) 



1 + + \ F £t 



is then 



(12) 



where is the relaxation rate of the nanomagnet at site i in the absence of the DDI field 
£ ddl and is given by the (intermediate-to-high damping) Neel-Brown expression 
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,(0) 



1/2 



TW7T 



(1 -h 2 ) (1 + h) e 



+ (l-h)e 



-a t (l-hf 



2a. 



1/2 



T~ s \/ 7T 



xT(ai,h), (13) 



with r s = (XjHa) . In the high-energy barrier approximation, cr^> l,/i = x/2er<Cl the 



function F,- reads 



22|] 



5 1 



(14) 



where A is the Landau-Lifshitz damping parameter. 

One should note that the relaxation rate in Eq. ( fl~3l) applies to the case of a magnetic 
field applied along the anisotropy easy axis, and can then be rigorously used only for a 
textured assembly, i.e., with all anisotropy axes parallel to the applied field. For an assembly 
with randomly distributed easy axes, one should use the (cumbersome) expression of the 



relaxation rate in an oblique field 13j, |23j . In the present calculation we ignore this effect 



and use expression ([13]) for all moments in the assembly and then average over the direction 
of the anisotropy axes [see discussion in section HVj . The same approximation was used in 



Ref. 



21 



while the calculations of Ref. [j of the ac susceptibility did not require a finite field 
and hence expression (fi~3l) was used at h = 0. 

The next step consists in substituting for £j in Eq. (fl2"l) the expression given by Eq. (|TT|) 
averaged over the spin and anisotropy orientations. Averaging over the spin orientations 
yields [see appendix lA] 

' '«ll) = ^ e «» 

(15) 

/A2 \ _ fej»i) 2 A 

where [see Eq. (IA7j) et seq.] 



A, 



6 + 3 ^ o 



L I.J 



Therefore, the relaxation rate of a weakly interacting particle containing n,i Bohr magne- 
tons, embedded in a polydisperse assembly, can be written as 



r, * if > [i 



(16) 
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where we have collected the DDI contributions in 

In the case of randomly distributed anisotropy easy axes one obtains [see appendix [B] 

2 — ^2n 2 : 

e . T> ■ V ■ e = •=> - 

ij j ij 

9 1 



and 

^ (£^) 2 1 + F { 



e*. (is) 



3 2 

Now that we have the expression for the relaxation rate that includes the DDI contribu- 
tion, we may study the effect of the latter on the critical volume V c , or the corresponding 
number n c of Bohr magnetons, which is defined by the equation [see Eq. (jSJ) et seq.] 



r(n c ), (19) 



where v m is the measuring frequency. The problem then is to determine how n c (or V c ) 
changes in the presence of DDI, recalling that it is a function of temperature, field, and 
other experimental conditions such as v m . For this we combine Eqs. (fl3| [T6l [T9|) to obtain 

y^^v^n*,h) [i+s]. 

Now, using Oi = a ni [see Eq. fl2])] and dropping the index i, we rewrite this equation as 

1 In n + In T(a , n, h) ~ In _ sfo, n). (20) 

2 V 2^/ao / 



Since this equation has been derived in the case of weak DDI we may seek its solution 
n c as an expansion in terms of the DDI coefficient ^ [see discussion of the validity of this 
perturbation in section IIV] . Indeed, inserting n c ~ n c0 + 5n c in (120]) and expanding around 
n C) o, which is the solution of Eq. (|20l) without the DDI term 5, we obtain the following 
expression for Sn c 

5r ^ = -- —, — rr =ri ( 21 ) 

2 - oo^on^o ft) + n Cj0 ^ 

where a c = aon c ^ and 

${(Tc,h) = — — . 

(pi{cr c ,h) 
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FIG. 1: (Color online) bn c versus the dimensionless parameter k (k = a/D m , a being the inter- 
particle distance and D m the mean diameter) for the damping parameter A = 0.01,0.1,0.25 and 
temperature T = 15 K. 



with 



<p n (a c , h) = (l + h) n e~ a ^ 1+h)2 + (1 - h) n e-^ {1 - h)2 

5 1 



-K,o) = @o ^~ 



"0 



and 



e 



26 f n c n 



1 - 
1 - 



8A 2 cr n C)0 
5 1 



16A 2 a n, 



c.O 



2 _ 2 -T-, 



(22) 



where 5 is the standard deviation of the volume distribution and n m = n e &2 the mean 
number of Bohr magnetons corresponding to the mean volume of the assembly; being 
the number of Bohr magnetons contained in the volume Vq [see Eq. Q]. 1Z ~ 16.8 jjj, [nj 
for a simple cubic lattice [see discussion in section UVj . 

The function $(cr c , h) decreases monotonically from 1 to when h varies from and 
tends to 1, and is nearly independent of a especially for large h. This implies that when 
the applied field increases 5n c increases (in absolute value) and thereby the effect of DDI is 
enhanced [see further discussion in the next section]. 
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FIG. 2: (Color online) ZFC magnetization of a non- interacting assembly for various field values. 
5 = 1. 

IV. RESULTS AND DISCUSSION 



In Fig. [T] we plot 5n c as a function of the dimensionless parameter k = a/D m with varying 
damping parameter A. We see that for the relatively small values of A, 5n c becomes negative 
and decreases with the increasing intensity of DDI (or decreasing k). This means that the 
critical volume, separating the dominating populations of blocked and superparamagnetic 
particles, decreases in the presence of DDI. This can also be seen from Eq. (|T^|) upon noting 
that the function F, and thereby the contribution of the transverse component of the DDI 
field in Eq. (II 2p . changes sign upon varying the damping parameter. A similar behavior was 



observed in Ref. 



24( (see Fig. 2 therein) where the energy barrier distribution for coupled 



Co particles was computed (using a different approach) as a function of their concentration. 
It was shown that the DDI induce an increment of the amount of small barriers, responsible 
for faster decay. In our case, this is equivalent to the decrease of V c under the effect of DDI. 



Obviously, all the curves (for different A) tend to zero for large k (absence of DDI). The 
limit between the negative and positive Sn c is given by the equation Ai imit = ^/5/(8a n Cj0 ), 
which implies that this limiting damping depends directly on temperature and on all other 
parameters via n c . 

Using the volume V c corresponding to n c ~ n Cj o + Sn c , with n c being the solution of 
Eq. (120]) without the DDI term H and 5n c given by fl2Tl) . we compute the ZFC magnetization 
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FIG. 3: (Color online) The temperature T max (H) as a function of the applied field for a) maghemite 
and b) cobalt particles for the same parameters D m and 5 as in Fig. [2j "Free" stands for the non- 
interacting assembly and the other curves are for the interacting case with mean interparticle 
distance a = kD m . The damping parameter is A = 0.01. 

according to Eq. fllOp . In Fig. [2] we plot the ZFC magnetization thus obtained as a function 
of temperature for various values of the applied field of a polydisperse assembly of non- 
interacting maghemite nanoparticles with mean diameter D m = 5 nm and standard deviation 
5=1.0 and random anisotropy. Apart from the obvious bell-like shape and the vertical shift 
of the maximum with the increasing field, we can see that the position of the maximum 
changes with the field in a non-monotonic way. 

In Fig. [3] we plot the position of the maximum of the curves in Fig. [2] as a function 
of the applied field, i.e., T max (if), for various values of the inter-particle (center-to-center) 
distance, for two substances. First of all, we observe that indeed the effect of DDI is to change 
T mayL (H) from a bell-like curve with a maximum to a monotonically decreasing function, and 



this compares well with the experimental results [see Fig. 1 of Ref. 



Ref. [si, the effect of DDI is not merely to change the potential energyscape, as was argued 



. As was stressed in 



in many previous publications 
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25 



26], but also to introduce a transverse field that 
induces saddle points in the potential jzj. This in turn makes the relaxation rate quite 
sensitive to the damping strength, and for relatively low damping, as is the case in Fig. [3] 
(A = 0.01), the probability of switching increases. In addition, if a magnetic field is added 
with increasing intensity, the energy barrier is lowered and the magnetic moments switch 
at lower temperatures. This concomitant effect makes the relaxation time shorter and for 
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a given observation time, the temperature T max (if) increases. This is of course compatible 
with the effect, discussed above, that DDI reduce the population of superparamagnetic 
particles in favor of the blocked ones and this leads to a larger T max (H). In summary, the 
critical volume V c (or equivalently n c ) given by Eq. ( 1201) in fact corresponds to the "critical" 
field at the maximum of T max (if) which, in the presence of the DDI field, separates: i) the 
low-field regime: the population of blocked particles becomes more and more dominant, 
which leads to the increase of the "average" blocking temperature, and ii) the high-field 
regime: the population of superparamagnetic particles takes over thus leading to a decrease 
of the average blocking temperature. Furthermore, increasing the intensity of the DDI 
field increases both its longitudinal and transverse components. While the longitudinal 
component contributes to the increase of T max the transverse component shifts its maximum 
towards low magnetic fields. 

Finally, we recall that the work of Ref. [§1 showed that, in the absence of DDI, the max- 
imum of T max (if) can be explained by the nonlinear variation of the superparamagnetic 
contribution to the ZFC magnetization with the applied field. In the presence of DDI we see 
that the contribution of superparamagnetic particles is reduced in favor of blocked particles 
and this leads to the disappearance of the maximum. 

The analytical expression (T2~T!) derived for 5n c is valid, in principle, for £d n m <C 1 [see 
notation in Eq. (j2J)], where n m is the mean number of Bohr magnetons in the assembly. This 
condition is equivalent to | 5n c /n Cj o |<C 1 [see Eq. (f2~Ti) ] which leads to the condition on k 



Then, using the physical parameters of Fig. [3] (left) with D m = 5 nm, n m ~ 2117, taking 
T = 20 K, H = 100 Oe and computing n Cj0 from Eq. fl2DD without the DDI term, Eq. (125)1 
yields k ^> 3, which is a reasonable condition. Accordingly, in Fig. [3] we took k > 4. 

Regarding the expression of the relaxation rate f[T3"j) a few remarks are in order. As 
we said earlier, rigorously, this expression applies to the case of a textured assembly with 
all easy axes parallel to the applied field. For random anisotropy, one should employ a 
numerical procedure for computing the relaxation rate with an oblique field, as is done 




(23) 



where 



o c = cr n Ci0 , f c = 1 



5 1 
16A^ 



$ c = $(cr c ,h). 
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in, e.g., Ref. [13|, [23|. However, we recall that the main objective of the present work was 
to i) understand the effect of DDI on the T max (if ) curve and ii) provide relatively simple 
(approximate) expressions including the DDI contribution. In addition, as the latter is only 
possible using perturbation theory which assumes weak DDI, a transverse magnetic field 
would be dominating and the subtle effect of the DDI field transverse component would not 
be easy to disentangle in a non ambiguous manner. On the other hand, we have shown 
here that the disappearance of the maximum of T max (H) is mainly due to the effect of the 
transverse component of the DDI field on the relaxation rate of the magnetic moments. One 
could then ask why the transverse component of the applied magnetic field does not play the 
same role in the case of a free assembly. The main reason is that the applied magnetic field 
has a static effect while the DDI provide a field that changes dynamically with temperature 
and other physical parameters related with the dynamics of the system. In this respect, we 
wish to make a connection with the work [^tJ where the effect of exchange interaction on 
the relaxation rate of a two-spin system was (semi-analytically) investigated by the kinetic 
Langer's theory. It was shown that, in the weak coupling regime, when the first spin starts 
its switching process and arrives at the saddle point, the orientation of the second spin 
undergoes some fluctuations creating a small transverse field that increases the switching 
probability or the relaxation rate. 

Furthermore, we would like to point out that in the present work we consider an as- 
sembly of nanoparticles placed at the sites of a regular (simple cubic) lattice, with varying 
inter-particle distances. The effect of changing the lattice structure or equivalently, the dis- 
tribution of the vectors [see Eq. flH]) et sea.], is to modify the lattice sum in Eq. ( 1221) . and 
thereby to change 5n c in Eq. (pTj) . In Ref. y the lattice sum 1Z (oc O , for a monodisperse 
assembly) takes the values 16.8, 14.5, 14.5 for simple cubic, bcc, and fee lattices, which leads 
to weaker DDI. On the other hand, in many realistic samples the position of the particles 
on the hosting matrix is random, and for a given concentration, the precise variation of the 
lattice sum Oo should depend on the particular form of the particles spatial distribution 
function, which may dramatically change in the presence of aggregates, chains, and the like. 

We also found that changing the volume distribution 5 width has only a quantitative 
effect on T max (H), very much similar to the results presented in Fig. 4 of Ref. [8j for a free 
assembly. Hence, in both interacting and non-interacting assemblies, we find that when the 
volume distribution becomes in magnitude and slightly flattens. 
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In addition, changing from maghemite to cobalt particles, which is mainly equivalent (in 
the present approach) to changing the anisotropy constant by an order of magnitude, has a 
quantitative effect on the curve T max (H, ^) curves but the qualitative features remain the 
same. 



V. CONCLUSION 



We have investigated the effect of (weak) dipolar interactions on the zero-field-cooled 
magnetization by computing their contribution to the longitudinal relaxation time. We have 
shown that the effect of the dipolar interactions is to lower the critical volume of the assembly 
which separates the dominating populations of blocked and superparamagnetic particles. 
More precisely, it is demonstrated that the maximum of T max (if) shifts towards low values 
of the applied field as the intensity of the dipolar interactions, or equivalently the sample 
concentration, increases. This result is in good qualitative agreement with experiments on 
both maghemite and cobalt nanoparticles. We finally emphasize the important role played 
by damping in the presence of a transverse field provided here by the dipolar interaction. 
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APPENDIX A: SPIN AVERAGES 

We compute the average of the square of the longitudinal and transverse components 
with respect to the local easy axis ej of the effective field Q [see Eq. ( fl"2l) ] comprising both 
the applied magnetic and the DDI fields with the condition \(\ <C 1. The magnetic field is 
included here only for completeness and is dropped in the final expressions obtained in these 
appendices. In fact, the calculation of the ZFC magnetization and thereby that of T max (if) 
requires the full range of this field, and for this reason in the calculation of T max (H) the 
magnetic field is included exactly in the relaxation rate [see Eq. ( Fl3l) ]. Nevertheless, the 
expressions of the longitudinal and transverse components of the effective field Q obtained 
here may be used in the range of small magnetic fields. 

The effective local field Q reads 

Q = Xie h + £ d 7ii Vij ■ Sj, (Al) 

3 

where Sj = rijSj. Its longitudinal component is then defined as 

Ct,|| = d ■ e i = Xi( e h ■ ei) + £ d n, ^2 e i ' ■ Sj. (A2) 

3 

and thus 



3 

+ {idnif 22 ' Vij ■ Sj) (ei ■ V ik ■ S fc )) 

3 k 

where the average () is defined with respect to the Gibbs probability distribution containing 
only the anisotropy term. 

Now, using the following formulas 

hi - s^s* + Sffetf 



(«r> = °. { s i4 



with 
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28] 



(2l+l)\\ \ 2 > T ■ ■ ■ , "i ^ ±, 



Sil((Ti) ~ < 



1 _ 1(1+1) , 



Oi > 1 



(A3) 



we obtain 



Kll)o = ^ ( eh ' e ^ 2 + 2x ^d.rii (e h ■ a) ^2 e i ' V ij ' ( S i)c 
+ (^i) 2 E E ' ^ ' S ^ ' ^* ' Sfc ^o 



j k 



x\ (e h ■ e,) 2 + n) [(1 - S j2 ) (et ■ V tJ ■ V VJ ■ e % ) + 3^ (e, ■ V t 



The transverse field is given by (Ci±) — (Ci)o ~ ' w ^ n 

<c 2 > = ^ 2 + (WE^<(^) 2 >o 



^2 n ) K 1 - ^2) (e* ■ % ■ % • e,) + 35^-2 (e* ■ V tj ■ e^) 2 ] 



(l - (e h ■ ei) 2 ) 



3 <(A i s i ) 2 > - (1 - S j2 ) (e, ■ V tj ■ Vij ■ *) 



-3S j2 (ei ■ ■ v.j ) 

Next, using Eq. (1A3jl . we compute the first term in the square brackets as 



:^,s,r) = ((Vi jS ,)(Vi jSj )) = Y, v if v iH s ^] 



a/3"/ 



a - « E « + 3 ^ E ( e ^J 

a/3 a/37 



Let us now compute these two terms. 



E« = E 



3egeg - i"" 3egeg - S"" 

X 



a/3 



a/3 

9 £4444 +i£^ 



r 



<> / j e ii e ii e ii e 



l a/3 



» a/3 



6 Z-r- W *J 



'J 



a: 



Next, 
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E 



ij ij 



944 E (44) - 3 4 E 4^ - 3 4 E 4^ + E rt 

a a a a 

- <h + 4-^ = 4 



y 



then 



a/3"/ 
1 

y 



7 



Recapitulating, we have 



!(% s j) 2 >( 



6 6 1 6 



r 



6 1 

l y y 



Therefore, inserting all results back into Eq. (1A5[) . we finally obtain 

(Cd^) 2 

(1 - ^2) (ej • Pjj • • ej) - 3S j2 (e$ • • e^- 



Gathering the results for both longitudinal and transverse components for the effective 
field, we write 



3 



(A7) 



[l - (e, • e,,) 2 ] + Y>f 



^7 + ^r.Sj'&ij - @y 



where we have introduced the notation 



9y = (1 - 5 i2 ) (e, • V iS ■ Vij ■ Bi) + 3S j2 n^ 
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APPENDIX B: AVERAGING OVER ANISOTROPY 



The general expressions for the longitudinal and transversal fields can be simplified in 
some relevant situations. For a textured assembly (with parallel anisotropy axes) we set all 
the e, parallel to e. For a system with randomly distributed anisotropy axes one replaces 
expressions involving /(e^) by integrals J d 2 e /(e) = /, and uses (e • vi)(e • v 2 ) = |vi ■ v 2 . 

In the present work we only deal with random anisotropy. In this case we have 

J d 2 e (e h ■ eif = i, 
so that for the longitudinal component we obtain 



3 



i 2 , 

" 3 * 3 ^ r% 
where we have used the averages 



1 — - 2 2 



6j ' T-^ij ' T-^ij g 3 [Qi • Gij) +1 ,, • 

[ e i ' ' e j] = ~6"3 (ej ■ eij) + 1 = - -Q-, 

ij ij 



Qij = (1 - S j2 ) (e, • V l3 ■ Vij ■ + ZS j2 n% 
2 „ , 2 1 „ 2 



-g- (1 - S j2 ) + 3--£-S j2 - 



'ij 'ij r ij 



For the transverse component, 



2 2 , {idUif 2 6 



3 r.6 

ij 



2 2 , V 4 ^ 2 



r 6 

ij 
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with 



1 

ij 



and thereby 



3 1 3 ^ r? 



Finally, for random anisotropy we have 



\ _ 1 _2 , fe" t ) 2 V- 2n i 
kll/ -3^+ 3 2_.^> 

772 \ _ 2 2 , {^nrf V 4n| _ / , 



(Bl) 
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